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Abstract 

A new exactly solvable alternative to the Calogero three-particle model is proposed. 
Sharing its confining long-range part, it contains the mere zero-range two-particle 
barriers. Their penetrability gives rise to a tunneling, tunable via their three inde- 
pendent strengths. Their variability can control the removal of the degeneracy of the 
energy levels in an innovative, non-perturbative manner. 
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1 Introduction 



In nuclear, atomic, molecular and statistical physics, the enormous popularity of the 
Calogero's solvable models [|l| reflects a fairly realistic form of their separate two- 
body interactions. They combine a long-range quadratic attraction with a strong 
short-range repulsion. In its full generality the model binds an arbitary number 
of particles moving along a straight line and its exact solvability is a consequence of 
certain deep symmetries of its partial differential Schrodinger equation 0. 

For the sake of brevity, let us pay attention just to the Calogero's first nontrivial 
three-body Hamiltonian 
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We may immediately see that from a purist's point of view, its two-body barriers 
are not penetrable and, in this sense, they do not admit any mutual exchange 
of particles. For this reason we proposed recently a solvable modification of this 
model where a tunneling has been rendered possible at an expense of a loss of the 
Hermiticity of the Hamiltonian We employed a complexification based 

on the shift of the coordinates x^'^^"'''^ — > x^^°™-P'-'^^\t) = t — ie{t^). This replaced the 
Hermiticity of the Hamiltonian by its mere commutativity with the product of parity 
V and complex conjugation T [Q. We have shown that the spectrum remained real 
in a way attributed to the VT symmetry in the related literature ||^. 

In the latter innovative few-body implementation of the idea of the VT sym- 



metrization, the Calogero's strongly singular real barriers with x 
were all replaced by the complex and e— dependent expressions 
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In the present short note we intend to describe and analyze an alternative scheme 
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which would admit a tunnehng. In essence, we shall start from the same leading-order 
formula (|^) but succeed in returning to a Hermitian Hamiltonian. 



2 The new model 
2.1 Inspiration 

Our new proposal has been inspired by formula (01) and by its Hermitization 
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e > 0. (3) 

e 



_x — le X + le 

The limit e — would then just reproduce the Calogero's model, one of the specific 
and most inspiring features of which lies in its separability in the three-body case. 
This is based on the re-parametrization of the coordinates Xi—X2 = \/2 p sin0 with 
p G (0, oo) and G (0, 27r) etc (cf. Figure 1). In the units 2m = h = 1 and for the 
equal and non- negative coupling constants (?! = (72 = S's = fl' > in eq. (|I]), such a 
change of variables reduces the Calogero's partial differential Schrodinger eqution to 
the mere ordinary Sturm - Liouvillean problem 

{-^, + Tr^frrrl ^ = ^"m, m = V^(2vr) = O. (4) 



2 sin^ 3 > ^ 

It solution generates the well known (viz., Laguerre times Jacobi) polynomial wave 
functions as well as the related equidistant spectrum with a gap, ~ + const, 



n = 0,2,3,4, 



2.2 Main idea 

Our attention is attracted by the tilded expression S^ix) in eq. (^. We shall try 
to re-read it as an approximate e ~ form of the well known Dirac delta function 
So{x) = 6{x). Thus, picking up, say, the first part of the two-body barrier in eq. (|I]) 
in its regularized form 
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we may insert the appropriate definition X2 — = \plp sin(0 — 47r/3) and get, at 
the very small e ^ 0, 



f\^e{x2 - 3:3) ^ ^ 
P 



In place of (^) it gives 



TT 
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27r 
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i.e., an approximative innovation of the Calogero's angular Schrodinger equation. 



2.3 Interpretation 

Due to the purely intuitive form of the above "derivation" of eq. (^, one has to be 
very careful in all questions related to its possible physical as well as mathematical 
interpretation. At the same time, the use of the delta-function-shaped potentials is 
quite common in practice as it makes many systems exactly solvable in the limit 
£ ^ 0. 

Let us start our further analysis of eq. by picking up some three constants 
Gk > 0. Then we re-define fk = fk{p) = V^Gk/ p and gk = £/fc(p)/7r. In the generic 
case with Xi —Xj 7^ 0, the limiting transition to e = makes all the three centrifugal- 
like forces vanish, gk 0. Such an observation is compatible with the philosophy of 
using just the contact barriers in our angular Schrodinger equation. 

In the second step we have to re-analyze the role of the overall singularity in the 
origin p = 0. In principle, we might admit and consider many different singularities 
there 0. Here we shall simplify our life by the most straightforward postulate that 
the point p = is just "regular" (i.e., "ignored" by our Hamiltonian) , and that 
all our radial wave functions are simply vanishing in the origin. Under such an 
assumption, the final p 7^ form of our new Hamiltonian if ("'^'") can be assigned a 
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formal representation 
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(xi - Xa)^ + {X2 - x-if + {x-i - X]f 
+ VLxb{x2 - X'^ + 1^25(a;3 - + fi3'^(a;i - ^2 

where the strength of the two-body contact terms has the three-body character and 
weakens with the distance of the detached spectator particle, 
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We may summarize that our Hamiltonian represents a new three-particle 

model which is exactly solvable. 



3 Solutions 



Our new angular Schrodinger equation 
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f)-(^-f) 



+ G3 [5 (0 - Tt) + 5 (0)] V^(0) = 4 V^fc(0) (7) 

describes the 0— dependent part of the three-body wave function and determines its 
energies via the Calogero-type formula, 

/3 



m,k 



■uj(2m+l + K,k), m,k — 0,1,... . 



(8) 



Within the open subintervals of 4> G T^j — {{j — 1] 7r/3, j 7r/3) we have a free motion 

^(0) = ^,(0) = Aj sm{K [0 - (j - 1) 7r/3]} + B, cos{k [0 - (j - 1) 7r/3]}. 

On the boundaries we have to match the separate local wave functions using the 
standard rules 



V-j ( ^ ) = Bj+i, 



jir 



j = l,2,...,6 
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with Pj{K) — (3j+3{K) — Gj/k and j — 1,2 or 3. The first five apphcations of these 
rules define all the wave functions V'i(0) iii terms of the initial one, say, V'i(0)- The 
sixth step becomes a selfconsistent matching which guarantees that the global wave 
function remains single- valued. In terms of the two auxiliary matrices 

R = R{k) 
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y — sin(fi;7r/3) cos(fi;7r/3) y y 
our matching conditions factorize in the two independent two-dimensional forms 



{u±i)\ Uo 




where U = U{k) = {I + l3zL)R {I + i32L)R {I + PiL)R. The pertaining two alternative 
two-dimensional secular equations 

det[W(fi;W) ±/] = (9) 

define all the roots k'^^ > with the sign ambiguity a — ±1 and the angular quantum 
number /c = 0, 1, ... in implicit manner. 



3.1 Toy model with the single barrier 

Let us fix d = G2 = and vary the strength G = G3. With (3 = G / k^"^ in the 
matrix 

U^''^\n) = (7 + /3L) • 



cos KTi sm «; TT 

— sin K TT cos K TT 



our determinantal secular equation reads 

G sin K TT 
cos kt: -\- a -\ = 0. 

2k, 

Ki a — —1 and any integer shift of — 2k -\- 8 with k — 0, 1, . . ., this secular 
equation gives us two roots 5(a), S{b) e [0, 1). The smaller one is trivial, 5{a) — 0. 
The second one is uniquely specified by the implicit formula 

7T6{b) G 



tan ■ 



4:k + 26{b)' 
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Also for cr = +1 and /t*-"*"^ = 2k + 1 + 5 we get 5{a) = and the similar relation 

7T6{b) _ G 

~ 4k + 2 + 25(6) ■ 

In the limit G ^ of the vanishing barrier we get the standard square-well solutions 
with the correct degeneracy 6{b) —>■ 6{a) at every k. 

The sign of the coupling depends on our choice, G G (—00,00). This provides 
an interesting counterpart to the Calogero model where the barrier cannot be too 
attractive [|T^. In the limit G — > 00 of the very strong repulsion, one returns to the 
Calogero-like case characterized by the impenetrability of the barriers. 

3.2 Three equal barriers 

When we choose Gi = G2 = Gs = G and abbreviate {I + j3 L)R = TZ{k) with 
/3 = G/k, our secular equation (^ may be factorized into the two separate conditions 
at both the quasi-parities a = ±1, 

det[7^(K) + al] = 0, det[n^{k) - an{k) + /] = 0. 

The first one parallels our preceding toy model and its solution is immediate, 

K{a) = 3N, K{b) = 3N + 35{b), N = 0,1,.... 

The implicit definition of the pertaining shifts A(6) G (0, 1) has almost the same form 
as above, 

7r6(b) G 
tan — — = — . 

2 6A^ + 66{b) 

In the tilded case, the identity det(7?.^ — aTZ + /) = (a — trT^)^ holds as long as 
detTZ = 1. This reduces the second secular condition to the trigonometric equation 

kn G k-K 
2 cos 1 sin — = cr. 

3 k 3 

Its solutions must be sought numerically, giving k = 1.367840720, 2.199769250, . . . 
at G = 1 (cf. Figure 2). 
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4 Summary 



One of the most striking features of the model of Calogero is the fully impenetrable 
character of its two-body repulsive barriers. They divide the phase space in the 
six independent subdomains (cf. Figure 1). In a way, this absolute impenetrability 
of the barriers is in its effect responsible for the exact solvability of the Calogero's 
model. 

This role of the Calogero's g barriers is partially weakened when they become 
attractive, i.e., > (7 > — 1/4. An extension of the Calogero model to this transition 
region has been discovered and described by Gangopadhyaya and Sukhatme In 
their construction the contact terms also appeared as a formal means of preservation 
of the solvability. As a consequence, the spectrum only consisted of the two shifted 
sets of equally spaced energy levels. Such a type of the modified Calogero's spectrum 
re-appeared also in our recent non-Hermitian construction 

In the present letter we were able to get rid of the Calogero's power-law 0{x^'^) 
barriers completely. Our alternative way of introduction of the contact terms enabled 
us to treat their couplings as independent parameters. Our key point is that these 
barriers are "thin", i.e., partially penetrable. This represents their main phenomeno- 
logical appeal. At all their finite (and, in fact, both repulsive as well as attractive) 
couplings Gj, their free variability might prove appealing in many phenomenological 
considerations. 

In contrast to the Calogero model characterized by the absolute absence of tunnel- 
ing, all our particles are permitted to jump over one another. In certain applications 
to few body systems of quark type, this could improve our intuitive insight and build 
some analogies with the motion in more dimensions. After all, the use of the har- 
monic two-body forces with an additional, contact "local spike" might also extend 
the advantages of the exact solvability quite easily beyond the traditional domains 
in the theoretical nuclear physics. 
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1 Introduction 



In nuclear, atomic, molecular and statistical physics, the enormous popularity of the 
Calogero's solvable models [1] reflects a fairly realistic form of their separate two- 
body interactions. They combine a long-range quadratic attraction with a strong 
short-range repulsion. In its full generality the model binds an arbitary number N 
of particles moving along a straight line and its exact solvability is a consequence of 
certain deep symmetries of its partial differential Schrodinger equation [2] . 

For the sake of brevity, let us pay attention just to the Calogero's first nontrivial 
three-body Hamiltonian 



2m 



g2 g2 g2 

+ 



dxi^ 8X2^ 9X3^ 
+ ^ [(Xi - X2f + {X2 - X^f + (^3 - Xif\^ + 
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(X2-X3)2 ' (X3-Xi)2 ' {xi - X2Y \ ' 

We may immediately see that from a purist's point of view, its two-body barriers 
are not penetrable [3] and, in this sense, they do not admit any mutual exchange 
of particles. For this reason we proposed recently a solvable modification of this 
model where a tunneling has been rendered possible at an expense of a loss of the 
Hermiticity of the Hamiltonian i^^^"') [4]. We employed a complexification based 
on the shift of the coordinates x^'"'^"^^ ^(complex) ^^-j _ ^ _ xhis replaced the 

Hermiticity of the Hamiltonian by its mere commutativity with the product of parity 
V and complex conjugation T [5]. We have shown that the spectrum remained real 
in a way attributed to the VT symmetry in the related literature [6]. 

In the latter innovative few-body implementation of the idea of the VT sym- 
metrization, the Calogero's strongly singular real barriers with 
were all replaced by the complex and e— dependent expressions 

{X - ief " X2 + £2 + (2;2 + ^2)2 + ^^l^ )" (" ) 

In the present short note we intend to describe and analyze an alternative scheme 
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Figure 1: The choice of coordinates for three particles 

which would admit a tunneling. In essence, we shall start from the same leading-order 
formula (2) but succeed in returning to a Hermitian Hamiltonian. 

2 The new model 



2.1 Inspiration 

Our new proposal has been inspired by formula (2) and by its Hermitization 



1 

2%^ 



X — IS X + IS 



TV 



Ss{x), 



£ > 0. 



(3) 
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The limit £ ^ would then just reproduce the Calogero's model, one of the specific 
and most inspiring features of which lies in its separability in the three-body case. 
This is based on the re-parametrization of the coordinates xi — X2 = \/2 p sin (j) with 
p e (0, oo) and (j) e (0, 27r) etc (cf. Figure 1). In the units 2m = h = 1 and for the 
equal and non-negative coupling constants fi'i = = = fl' > in eq. (1), such a 
change of variables reduces the Calogero's partial differential Schrodinger eqution to 
the mere ordinary Sturm - Liouvillean problem 
f (f 9g ' 



V'(O) = ^{27^) = 0. 



(4) 



V dcjy^ 2 sin^ 3 

It solution generates the well known (viz., Laguerre times Jacobi) polynomial wave 
functions as well as the related equidistant spectrum with a gap. En ^ n + const, 
n = 0,2,3,4,.... 



2.2 Main idea 

Our attention is attracted by the tilded expression Se{x) in eq. (3). We shall try 
to re-read it as an approximate e form of the well known Dirac delta function 
So{x) = S{x). Thus, picking up, say, the first part of the two-body barrier in eq. (1) 
in its regularized form 

9i n c r \ r T^gi 



{x2 - x^y + 



flS,{x2 - Xs), 



fl 



e > 0, 



we may insert the appropriate definition X2 — Xs = y/2 p sin((/) — 47r/3) and get, at 
the very small e — )• 0, 



flS,{x2 - X3) R 

In place of (4) it gives 



47r\ 



fiP 



-G2 



i.e., an approximative innovation of the Calogero's angular Schrodinger equation 



(5) 
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2.3 Interpretation 

Due to the purely intuitive form of the above "derivation" of eq. (5), one has to be 
very careful in all questions related to its possible physical as well as mathematical 
interpretation. At the same time, the use of the delta-function-shaped potentials is 
quite common in practice [7] as it makes many systems exactly solvable in the limit 
e ^ 0. 

Let us start our further analysis of eq. (5) by picking up some three constants 
Gk > 0. Then we re-define fk = fk{p) = V^Gk/p and gt = e fk{p)/'^- In the generic 
case with Xi — Xj ^ 0, the limiting transition to e = makes all the three centrifugal- 
like forces vanish, — > 0. Such an observation is compatible with the philosophy of 
using just the contact barriers in our angular Schrodinger equation. 

In the second step we have to re-analyze the role of the overall singularity in the 
origin p = 0. In principle, we might admit and consider many different singularities 
there [9]. Here we shall simplify our life by the most straightforward postulate that 
the point p = is just "regular" (i.e., "ignored" by our Hamiltonian), and that 
all our radial wave functions are simply vanishing in the origin. Under such an 
assumption, the final p ^ form of our new Hamiltonian can be assigned a 

formal representation 



h 

2m 



Q2 Q1 Q2 



+ 



dxi^ dx2^ dxs"^ 

+ ^ (J^ [(Xi - X2f + {X2 - XsY + {Xs - Xif^ + 

+f^i^(a:2 - xs) + fl2S{x3 - Xi) + fl3S{xi - X2) (6) 

where the strength of the two-body contact terms has the three-body character and 
weakens with the distance of the detached spectator particle, 

"1 — 1 T; ^'2 — 1 r, "3 



Fl - X2\ \X2 - X3I |X3 - Xi\ 

We may summarize that our Hamiltonian represents a new three-particle 

model which is exactly solvable. 
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3 Solutions 



Our new angular Schrodinger equation 

^2 



+G2 



+G3 [5 (</. - tt) + 5 (</.)] ^(</.) = M^) (7) 

describes the ^—dependent part of the three-body wave function and determines its 
energies via the Calogero-type formula, 

^3 



E = Em,k = \l^(^{'^rn + l + Kk), m,k = 0,l,... 



(8) 



Within the open subintervals of G Vj = {[j — 1] 7r/3, j7r/3) we have a free motion 

i;{<P) = ^j{<t>) = Aj sm{K [cp - {j - 1) 7r/3]} + cos{k [0 - 0' - 1) n/3]}. 

On the boundaries we have to match the separate local wave functions using the 
standard rules 



B 



with (5j{K) = /3j+3(K) = Gj/K and j = 1,2 or 3. The first five applications of these 
rules define all the wave functions tpj{(t)) in terms of the initial one, say, tpi{(t))- The 
sixth step becomes a selfconsistent matching which guarantees that the global wave 
function remains single- valued. In terms of the two auxiliary matrices 

/ 



K[Aj+i + (3j{K)Bj+i], 



j = l,2,...,6 



R = R{k) 



cos(K7r/3) sin(K7r/3) 
^ — sin(K7r/3) cos(«:7r/3) j 













L = 






1 




.1 


J 



our matching conditions factorize in the two independent two-dimensional forms 

{'U±I) ^ =0 
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where U = U{k) = {I + l3sL)R {I + I32L)R {I + l3iL)R. The pertaining two alternative 
two-dimensional secular equations 

det[Z^(KW) ±/] = (9) 

define all the roots k!)^^ > with the sign ambiguity cr = ±1 and the angular quantum 
number A; = 0, 1, ... in implicit manner. 

3.1 Toy model with the single barrier 

Let us fix Gi = G2 = and vary the strength G = G3. With /3 = G/k^'^^ in the 
matrix 

U^'^y^K) = (7 + /5L) ■ 



^ cos K TT sin KTT ^ 



y — sm K TT COS ktt J 



our determinantal secular equation reads 



G SmKTT 

cos K TT + cr H = 0. 

2k 

At cr = — 1 and any integer shift of k'^ '' = 2k 6 with k = 0, 1, . . ., this secular 
equation gives us two roots 5(a), 5(6) e [0,1). The smaller one is trivial, 5[a) = 0. 
The second one is uniquely specified by the implicit formula 

7rS(b) G 

tan = — — — . 

2 4k + 25{b) 

Also for cr = +1 and k^~^^ = 2A; + 1 + 5 we get S{a) = and the similar relation 

7rS(b) G 

tan = — — . 

2 4k + 2 + 26{b) 

In the limit G ^ of the vanishing barrier we get the standard square- well solutions 
with the correct degeneracy S{b) — > 5{a) at every k. 

The sign of the coupling depends on our choice, G e (—00,00). This provides 
an interesting counterpart to the Calogero model where the barrier cannot be too 
attractive [10]. In the limit G — ^ 00 of the very strong repulsion, one returns to the 
Calogero-like case characterized by the impenetrability of the barriers. 
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Figure 2: Graphical determination of the roots k = \fE at Gi = G2 = = 1 

3.2 Three equal barriers 

When we choose G\ = = G->, = G and abbreviate (/ + (3 L)R = TZ{k) with 
(5 = G/ K, our secular equation (9) may be factorized into the two separate conditions 
at both the quasi-parities cr = ±1, 

det[7^(/^) + al] = 0, det[n^{k) - aTl{k) + 7] = 0. 

The first one parallels our preceding toy model and its solution is immediate, 

K{a) = 2,N, K{b) =3N + 3S{b), N = 0,1,.... 

The implicit definition of the pertaining shifts A(5) G (0, 1) has almost the same form 
as above, 

7rS(b) G 
tan = — — . 

2 6N + 6S{b) 

In the tilded case, the identity det{7l^ — aTZ + /) = (a — trT^)^ holds as long as 
deiTZ = 1. This reduces the second secular condition to the trigonometric equation 

kn G kn 
2 cos h — sm = a. 

3 k 3 

Its solutions must be sought numerically, giving k = 1.367840720, 2.199769250, . . . 
at G = 1 (cf. Figure 2). 
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4 Summary 



One of the most striking features of the model of Calogero is the fully impenetrable 
character of its two-body repulsive barriers. They divide the phase space in the 
six independent subdomains (cf. Figure 1). In a way, this absolute impenetrability 
of the barriers is in its effect responsible for the exact solvability of the Calogero's 
model. 

This role of the Calogero's gx'"^ barriers is partially weakened when they become 
attractive, i.e., 0>(7>— 1/4. An extension of the Calogero model to this transition 
region has been discovered and described by Gangopadhyaya and Sukhatme [8]. In 
their construction the contact terms also appeared as a formal means of preservation 
of the solvability. As a consequence, the spectrum only consisted of the two shifted 
sets of equally spaced energy levels. Such a type of the modified Calogero's spectrum 
re-appeared also in our recent non-Hermitian construction [4]. 

In the present letter we were able to get rid of the Calogero's power-law 0{x~'^) 
barriers completely. Our alternative way of introduction of the contact terms enabled 
us to treat their couplings as independent parameters. Our key point is that these 
barriers are "thin", i.e., partially penetrable. This represents their main phenomeno- 
logical appeal. At all their finite (and, in fact, both repulsive as well as attractive) 
couplings Gj, their free variability might prove appealing in many phenomenological 
considerations. 

In contrast to the Calogero model characterized by the absolute absence of tunnel- 
ing, all our particles are permitted to jump over one another. In certain applications 
to few body systems of quark type, this could improve our intuitive insight and build 
some analogies with the motion in more dimensions. After all, the use of the har- 
monic two-body forces with an additional, contact "local spike" might also extend 
the advantages of the exact solvability quite easily beyond the traditional domains 
in the theoretical nuclear physics. 
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Figure captions 

Figure 1. The choice of coordinates for three particles 

Figure 2. Graphical determination of the roots k, = y/E at 

Gi = G2 = (j3 = 1 
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